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A theory capable of explaining intrinsic and extrinsic tunnelling conductance in underdoped 
cuprates has been devised that accounts for the existence of two energy scales, their temperature 
and doping dependencies. The asymmetry and inhomogeneity seen in extrinsic (normal metal - su- 
perconductor (NS)) tunnelling and the normal-state gapped intrinsic (SS) conductance is explained, 
»' . as well as the superconducting gap and normal state pseudogap and the temperature dependence 

^ ■ of the full gap. 
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^ ' Since the discovery of high- Tc superconductivity in 1986 by Bednorz and Miiller fl*!, there has been a huge theoretical 

, effort to understand the mechanism behind it. A lanthanum barium copper oxide was the first compound displaying 
^ this phenomenon fH, now we know there are many compounds displaying high temperature superconductivity that 
^ involve copper and oxygen, these make up the cuprate family. Cuprates are distinguishable from conventional metal- 

lie superconductors by originating from the doping of the parent charge-transfer insulators. The superconducting 
^ parts are weakly coupled two dimensional doped layers held together by the parent lattice. Cuprates have unique 

1 t properties, as well as their high-Tc they also have two energy scales, or gaps: The BCS-like "superconducting" gap 
Ci (SG) present in cuprates and other related compounds develops below the superconducting critical temperature and 
^ can be seen by extrinsic and intrinsic tunnelling experiments as well as high-resolution angle-resolved photoemission 

(ARPES) experiments; there also exists another energy gap, the "pseudogap" (PC) which is a large anomalous gap 
that exists well above Tc- The PG phenomena was first observed in spin responses 0, @l ^^.d using scanning tunnelling 
spectroscopy (3) in underdoped YBa2Cu307_5. Shortly after, the same gap was observed through infrared measure- 
ments many experiments have since exhibited this PG. The first explanation of the gap was offered in the form 
of real-space preformed hole pairs [6| called small bipolarons which are bound together by a strong electron-phonon 
interaction (EPl). Since then many theoretical explanations have been proposed for the origin of the PG which can 
roughly be divided into two groups. The first of these groups argues that the PG originates from some order, either 
Qv^ ■ static or fluctuating. The second understands the PG is the precursor of the SG, and reflects pair fluctuations above 
p\J , Tc- Some of the theories from the first group see the superconducting state of cuprates as being the result of a 
■rj" doped Mott-insulator (for example Q, see Q for a review). In his resonating valence bond (RVB) theory, Anderson 
[ focuses on the ground state and low lying excitations, the origin of the PG is seen as the spin gap associated with 
T^lj- . the breaking of RVB singlets @ . It has been suggested that adding impurities to (or doping) cuprates could weaken 
' the order parameter (for example this order parameter could be antiferromagnetic spin fluctuations (loj ) and thus be 
, the cause of the PG. Some believe the PG could be the result of SU2 rotations (for example [lH), in the underdo ped 
• region which connect fluctuations of staggered flux states and d-wave superconductivity. Chakravarty et al (2001) [IJI 
^ , proposed a static orbital current state called a c?-density wave was the origin of the PG based on phenomenological 
grounds. It has been argued that the PG is a consequence of a spin density wave (SDW) or charge density wave 
(CDW) state [ll], or an interplay between the two [l3]- It has also been suggested that the PG could be the result of 
inhomogeneous charge distributions containing hole- rich and hole-poor domains [iH[3j or the cause of the SG and 
PG could be the inter-band pairing of an itinerant band and defect states [l3|- 

The second group bases the understandiiig of the PG and high-Tc superconductivity on pairing interactions, pre- 
formed Cooper pairs have been suggested [18| | where the pairing is not in real-space but instead in momentum-space. 
It has however been implied that the short coherence length of cuprate superconductors suggests they lie somewhere 
between the BCS limit of very large momentum-space pairs and the opposite case of small real-space pairs under- 
going a Bose-Einstein condensation (BEG) JJi]. The BCS-BEC crossover has been studied in detail, for example 
in Ref . a superfluid state is approached in a system of localised bosons (tightly bound electron pairs) in contact 
with a reservoir of itinerant fermions (electrons), it is assumed the spontaneous decay and recombination between 
the two species causes superconductivity and the PG is a consequence of this, opening up in the fermionic density 
of states (DOS). Attractive Hubbard models have been considered as the origin of superconductivity and the PG. 
For example, studies of the normal-state of the two dimensional attractive Hubbard model have been carried out 
using quantum Monte Carlo (QMC) calculations [2l[, also the excited and ground state properties of the two dimen- 



2 



k-integrated over 1st BZ 




La, SrCu04 



0.3 0.2 0.1 0.0 0. 

E-E, (eV) 



FIG. 1: (Colour online) Momentum integrated photoemission over the first Brillouin zone of La2-a:Sra;Cu04 [22|, showing no 
signs of the van Hove singularity. 



sional attractive Hubbard model have been studied using the conserving, self consistent T-matrix formalism in the 
intermediate coupling regime and at low electron concentration ^23i] . Other approaches emphasize the weak phase 
stiffness in underdoped cuprates which is a result of low superfluid density or superconducting carrier density and 
it leads to a suppression of Tc by phase fluctuations, this means the underdoped cuprates are characterised by a 
relatively small phase stiffness and poor screening [l8l |. A diagrammatic theory of the one-band Hubbard model was 
proposed, where the temperature of the onset of the PG is related to the scattering rate [24|. The effects of classical 
phase transitions on the quasiparticle spectra were contemplated in underdoped cuprates in the PG regime above 
Tc by taking into account mean-field d-wave quasiparticles that are semiclassically coupled to supercurrents induced 
by fluctuating unbound vortex-antivortex pairs 25f|. Another idea with incoherent d-wave quasiparticles suggests 
that when the phase-coherence length exceeds the Cooper pair size, a PG appears [2^, the phase fluctuations of a 
dj.2_y2 pairing gap in a two dimensional BCS-like Hamiltonian approach is thought to be the origin of the PG [l^l- A 
phenomenological theory was produced that allowed the modelling of the effect of local superconducting correlations 
and long-range phase fluctuations on the spectral properties of high-temperature superconductors by reasoning that 
the PG is connected to the character of the excitations that are responsible for destroying superconductivity 12811 . 

Femtosecond spectroscopy is a tool for studying the temperature dependence of gaps, for example in Ref.[29[ the 
temperature independence of the PG and dependence of the SG in Yi_a;CaBa2Cu307-5 and HgBa2Ca2Cu308+5 was 
found. Raman spectroscopies have also been able to find two energy scales [s^l- ARPES has provided valuable infor- 
mation about cuprates, ARPES performed on Bi22I2 UMI, Bi220I [ssj, LSCO (36|], LBCO [37|] and CaNaCuOCl 
[sst has verified the presence of two energy gaps in cuprates. 

Scanning tunnelling microscopy (STM) offers a powerful technique to look at the doping, temperature and spatial 
dependence of the DOS with high resolution. It is sensitive to the DOS near the Fermi energy and to a gap in 
the quasiparticle excitation spectrum. Extrinsic tunnelling experiments have left us with many questions regarding 
the properties of cuprates. STM tunnelling spectra exhibit an SG and PG [39l - l4l| whose origin currently remains 
unaccounted for, despite many ideas partially discussed above. STM results on single crystals of Bi22I2 (for example 
^3-01) and LSCO [isl, have demonstrated the temperature, doping and spatial dependence of the SG and PG. 
In particular, in NS tunnelling, Kato et al [46^1 found that the PG is not uniform in real-space and its spatial average 
increases with decreasing hole concentration in spite of suppression of critical temperature. On the other hand a 
smaller gap (presumable SG) is uniform across the sample and is less doping dependent. 

Intrinsic (superconductor-superconductor, SS) tunnelling experiments on small Bi22I2 [47l-[50| and LSCO 51 1 mesas 
have found sharp quasiparticle peaks at the SG and broad humps representing the PG [47|. The PG exists above and 
below Tc and can persist up to room temperature (47j . The advantages of intrinsic tunnelling are that it is a direct 
spectroscopic technique that avoids problems like surface deterioration [i^ , it probes the bulk electronic properties of 
samples, it offers high resolution whilst being mechanically stable so it is perfectly suited for temperature dependent 
studies of high- Tc superconductors [5Q|. Break junction experiments also exhibit the SG [H,!!!] in underdoped Bi2212 
samples, the coexistence of the SG and the PG is seen in Bi2201 [54] and Bi2212 [55| and quasiparticle energy gaps 
are found in Bi22I2 tunnelling spectra [56l. [57|. One puzzling observation in cupratcs is the significant asymmetry 



3 




FIG. 2; (Colour online) Ratio of the negative bias NS tunnelling conductance to the positive bias R = /jvs( — 100)//]vs(100), 
integrated from to =FlOOmeV respectively, carried out for some cuprate superconductors [H, [H^ over a wide range of atomic 
hole density. The two curves express the asymmetry you would expected to see from a Mott insulator (solid blue) and 
conventional semiconductor (dashed red) without electron hopping. 



between the negative and positive bias conductance in NS tunnelling, which means the direction of the tunnelling 
carriers affects the tunnelling conductance. One possible explanation of the asymmetry was offered by the van Hove 
singularity of the DOS, a singularity in the DOS would give a visible hump in the ARPES data, however this does 
not seem possible since this singularity is not present in the momentum integrated photoemission ^22], see Fig[T] 
and Ref.[60|. In the tunnelling spectra of conventional semiconductors or Mott-insulators asymmetry is expected. 
Consider a semiconductor where the number of electrons is twice the number of ions as each ion can accommodate a 
spin up and spin down electron. Removing X electrons from the sample leaves the number of electrons as 2N — X . 
When positive bias is applied to the sample, the electrons tunnel from the tip to the sample, application of negative 
bias gives tunnelling in the opposite direction. The probability of an electron tunnelling is, for negative bias (sample 
to tip), proportional to the number of electrons available this is 2N — X, and for positive bias (tip to sample) the 
probability of tunnelling is proportional to the number of holes available in the sample for the electrons in the tip 
to tunnel to, this is X. The ratio of the integrated negative and positive conductance is given by i? = (2 — x)/x, 
where x = X/N. Similarly we can consider a Mott insulator where the Coulomb repulsion is so strong that the 
number of electrons is equal to the number of ions as each ion can accommodate just one electron. Removing X 
electrons from the sample leaves the number of electrons as N — X. Following the same idea as for the semiconductor, 
we have i? = (1 — x)/2x (60| . This is a very basic formulation ignoring any electron hopping, it is used to give us 
an idea of the magnitude of the asymmetry we can expect to see in a semiconductor or Mott insulator. It can be 
seen in Fig|5]that although both insulators exhibit asymmetry, they do not account for the magnitude of asymmetry 
seen in STM experiments with cuprate superconductors without the consideration of disorder and matrix elements 
[6l| . Despite intensive research, a detailed microscopic theory capable of describing unusual ARPES and tunnelling 
data have remained elusive and so the relationship between the SG and PC has remained unknown. A detailed and 
consistent interpretation of SG and PC could shed light on the key pairing interaction in cuprate superconductors. 

High values of Tc and small isotope effect on Tc in optimally doped YBa2Cu306.9 led some authors to conclude that 
the pairing interaction between electrons cannot be mediated by phonons. However experiments [6^ [63| showed that 
a partial substitution of Yttrium by Praseodymium, or of Barium by Lanthanum leads to the isotope effect simulta- 
neously with the decrease of Tc'. either these substituted compounds have a different mechanism of superconductivity 
or the mechanism is always phonons and the absence of the isotope effect in YBCO is due to something else. In favour 
of the last option are the tunnelling spectra at higher voltages of NCCO [63| and BSCCO [6^ [6^. Evidence from 
the doping dependent oxygen effect (OIE) on Tc and the substantial OIE on the carrier mass suggests a strong EPI 
in cuprate superconductors, where lattice vibrations play a significant but unconventional role in high temperature 
superconductivity, see Ref.[67[ and references therein. 

In this paper we develop a theory of NS and SS tunnelling in the bosonic and cuprate superconductors extending 
our previous brief report [6l| . The theory is based on the assumption that the EPI is strong enough in cuprates and 
similar ionic charge-transfer insulators to form small mobile bipolarons, which has been convincingly supported by 
a number of experimental observations (68j . To clarify the terminology we first introduce bosonic superconductivity 
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(Section |TI]) . The results on the NS and SS tunnehing are presented in Sections IIIII and ITVl respectively. 



II. BOSONIC (BIPOLARONIC) SUPERCONDUCTIVITY AND CUPRATE BAND STRUCTURES 

The BCS theory is capable of successfully describing the superconducting properties of elemental superconduc- 
tors with a small EPI strength. The theory was modified in 1960 to give a strong coupling theory fld\ describing the 
properties of intermediate-coupling superconductors (the difference between a weak and strong-coupling superconduc- 
tor is given by the electron-phonon coupling constant, A [zi|). The mean-field BCS-Eliashberg theory is applied when 
the electron correlation length is large compared to the distance between them. The mechanism behind supercon- 
ductivity is the momentum-space pairing of electrons through electron-phonon interactions. It was first realised by 
Frohlich in 1950 that electrons could be attracted to one another through their interactions with phonons; he suggested 
superconductivity was instigated by EPI. The isotope effect observed experimentally verified Frohlichs proposition 
that EPI causes superconductivity. 

When the coupling constant is increased above A ~ 1, the kinetic energy of electrons becomes small compared with 
the potential energy from the local lattice deformation, thus all electrons in the Bloch band become dressed with 
phonons (for a recent review see [Z^l)- The electron becomes a quasiparticle, a small polaron, which can propagate 
through the lattice in a narrow (polaronic) band together with the lattice deformation. For a further extension of the 
BCS theory towards a strong interaction between electrons and ion vibrations, A > 1, it was predicted that instead of 
Cooper pairs, a charged Bose gas of tightly bound small bipolarons would be evident [t^ with a polaronic BCS-like 
high- Tc superconductivity in the crossover region [73 |. These bipolarons are real-space pairs of two electrons with 
their phonon cloud. 

Different from Cooper pairs in the momentum-space, the ground-state of the strongly coupled electrons and phonons 
is the real-space pairing of these single polarons into bosonic bipolarons where they form a condensate which can be 
described as a charged Bose-liquid on a lattice if the carrier density is small enough to avoid their overlap [gI]. In 
the superconducting state, if the temperature is finite, not all the polarons will condense and those that have not 
condensed interact with the _condensate through the same potential that binds them together. The single-particle 
Hamiltonian is described as /7S : 



(1) 



where = E^, — /z, Ey is the normal-state single polaron energy spectrum in the crystal field and disorder potentials 
renormalized by EPI and spin fluctuations and Ac^ — — Acp is the coherent potential proportional to the square root 
of the condensate density, Ac oc ^nc{T). The operators pi and p\, create a polaron in the single particle quantum 
state v and in the time reversed state u respectively. 

As in the BCS case, the single-particle energy spectrum is found by applying the Bogoliubov transformation to 
diagonalise the Hamiltonian, which is thus written: 



(2) 



where py = UyUy + w^/Jj, py = UyjSy " V y , ty = ^^^^ + A^^^ wlth u^, = i ^1 ± . This spectrum is different 

to the BCS quasiparticles because the chemical potential (/i), is negative with respect to the bottom of the single- 
particle band, /i = — Ap. A single-particle gap A, is defined as the minimum of the single-quasiparticle energy 



spectrum. Without disorder, for a point-like pairing potential with the s-wave coherent gap, Acfc ~ Ac, one has [75 



A(T) 



A2 + Ac(r)2. 



(3) 



The full gap varies with temperature from A(0) = y A^ + Ac(O)-^ at zero temperature to the temperature independent 

A = Ap above Tc, which qualitatively describes some earlier and more recent j5Q] observations including Andreev 
reflection in cuprates ((7S] and references therein). 

In this paper we adopt the parent band structure based upon a Mott or any semiconductor insulator structure. 
The copper band is split into two; an upper and lower band which makes the Mott insulator, the gap between the 
two is approximately 5 — %eV . The semiconductor oxygen p-band lies within this gap and its charge transfer gap 
is approximately 1 — 2eV ^ so we consider a parent lattice that is half Mott-insulator, half semiconductor. We have 
assumed the local density approximation and generalised tight binding ("LDA+GTB") band structure [tJ, see Fig|31 
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FIG. 3: (Colour online) Tlie"LDA+GTB" band structure is shown on the left with the impurity levels indicated, these impurity 
levels cause the bandtail in the cuprate superconductor DOS on the right. 



Here, following Ref.[78| we amend this structure with the impurity bandtails. When the cuprate is doped an 
impurity ion locally introduces a distinct energy level within the charge transfer gap. The random spatial distribution 
of impurities when doping causes a bandtail effect of the DOS, similar to that of a heavily doped semiconductor. 
This band structure explains the charge transfer gap, sharp quasiparticle peaks near (§. f ) of the Brillouin zone 



and a high energy waterfall observed by ARPES in underdoped cuprate superconductors [78j. Only the impurity 
states with binding energy below ^ — contribute at T = OK. It has been suggested that the band structure should 
be metallic due to Fermi arcs seen by ARPES which could form part of a large Fermi surface; however this does 
not take into account strong correlations and the existence of the charge-transfer gap over a wide range of dopings 
[ttIItqI. For tunnelling to be possible, a state must be occupied by a carrier and at the same energy level on the other 
material a state must be vacant. The probability of a single-particle state being occupied is given by the Fermi-Dirac 
distribution. We consider single-particle tunnelling only, this is in key with STM results as they measure one particle 
tunnelling only (rather than Josephson tunnelling where the tunnelling of pairs can occur). The tunnelling process 
is described by the standard perturbation theory, where the tunnelling Hamiltonians are perturbations, then the 
Fermi-Dirac golden rule (FDGR) is applied. 



To find the NS tunnelling Hamiltonian each different tunnelling scenario needs to be considered, FigUl Suppose 
on the left we have the metallic tip with undressed carriers as opposed to the polarons and bipolarons on the right, 
superconducting side. This means for tunnelling left to right we have the annihilation of a free carrier on the left 
accompanied by the creation of a polaron on the right. Alternatively we might have the annihilation of a carrier on 
the left and the annihilation of a polaron on the right with the creation of a bipolaron on the right. This can be 
expressed in terms of the Hamiltonian [76j : ("h.c." is the Hermitian conjugate, describing the tunnelling in opposing 
direction) 



Here and fej^, are the annihilation of a carrier in the metallic tip in state ly and the creation of a composed boson in 
the superconductor in state rj' respectively, N is the number of lattice cells. P and B arc tunnelling matrix elements 
respectively with and without the involvement of a bipolaron. Generally B > P, because the presence of an additional 
hole lowers the tunnelling barrier for an injection of the electron |76,]. Using the Bogoliubov coefficients 




III. NS TUNNELLING 
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FIG. 4: (Colour online) Cartoon demonstrating the two possible single-particle tunnelling scenarios. The first is the annihilation 
of an electron in the superconductor on the left and the creation of a polaron on the right, as described in the first part of the 
tunnelling Hamiltonian. The second illustrates the tunnelling process with the involvement of a bipolaron. For normal-metal 
to superconductor tunnelling, this is the annihilation of an electron in the metal and the annihilation of a polaron in the 
superconductor with the creation of a composed boson, this is described in the second term of Eq.Q. Energy is conserved in 
the tunnelling process. 



to replace the polaron operators with linear combinations of the quasiparticle operators yields: 

B 



Hns = {uy,al, + v^iPy^ Cy + -j= ^ h]^, {uy.^i,: - Vy-oi\)^ Cy -f- h.c. (6) 

Using the Fermi-Dirac Golden Rule 



W 1] 



W = ^ |({n} liJtunI {0})|'<5(i?„ - i?o), (7) 



where the initial state is and the final is n, yields 

Til 

H 



?2 



^ (1 + n^,) [ul,fy,FJ{E^, -^.-eV- e,0 + vl,{l - fyW{E,y - - eV + e,,)] , (8) 



N 



VF^s - E ["''/-'(I - + - ^-') + ^''(1 - - Fy)S{^y + eV + e,0] 



2^ 

N 



J2 [ul,il - Mil - F,)5{E,y - - el/ - 6,0 + - F,)5{E,y - - eV + e,.)] .(9) 



Here W™ and are transition rates in and out of the superconductor, fyi = l/{e'^"'/^'^'^ + 1) is the single 

quasiparticle distribution function, n^/ is the bipolaron (Bose) distribution function, F^, = 1/(6^"/^^"^ + 1) describes 
the distribution of carriers in the normal metal, V is the voltage drop across the junction and the bipolaron chemical 
potential in the left superconductor differs from the right one by 2eV . To find the current we use the equation 



which gives jGlj: 

Ins{V) 



h 



J2 K {F. - U')5{iu + eV-e,,) + (F, + /,y - 1) + eV + e,,)] 



(10) 



+ — ^ {ul, [F,U. - {x/2){l - - U,)] SiC. +eV + e,, 



where x/2 is the atomic density of composed bosons in the superconductor. The boson energy dispersion is neglected 
here for more transparency, assuming that they are sufficiently heavy for their bandwidth to be relatively small. Using 
~^ /^/'^^(O^C: X^jy' ~^ PJv(C')'^?'j ^-iid neglecting the energy dependence of the metallic DOS pm{0 since 
near the Fermi energy it is approximately a constant, we obtain: 



+ 



h 

2-KeB'^PM 



J d^PNia {uHn [F{^ - eV) - fie')] e') + v^O m - e^^) + /(^O - 1] 5{i + e')} 

d^ JdCpNia {uHO [Fi^ - eV)f{e') - | (1 - - eV) - /(e'))] 5 + e') 
+ v\a [fH - eV) (1 - f{e')) + I {F{i - eV) - /(e'))] + e')], 



(12) 



where pn{(,) is the normal-state single-particle DOS in the doped charge-transfer insulator. At zero temperature, the 
Fermi-Dirac distribution becomes a step function, F(^ — eV) — > Q{eV — ^) and /(e') = 0. To find the conductance, 
the current is differentiated with respect to the voltage and we have: 



<^NS 



2TreP^PM 

h 

2'KeB'^PM 



di J d^PNiC) [n'imeV - - e') + v^aSieV - ^SiC + e')] 



(13) 



Using the Bogoliubov coefficients defined earlier, Eq.([5]), we find: 



a^S oc e(et/-A,){ps(eT/)(2 + |) [pNiVieV)^ - ^l) + PN{-^{eVY - A^) 

'pN{~^{eVY~/^l) - pNW{eVY-/^l)\ } 



+Q{-eV - A,) [psieV) (l + |) [p^WieVf-^D + pn{-^ [eVf ~ Al) ^ 

+ (l - f ) [pN{-^{eVY-Al) - PNWieVY-^l)] } 



(14) 



where 0(x) is the Heaviside step- function, ps{E) — E / ^ E"^ — A^ for s-wave symmetry of the coherent gap which does 
not depend on the quantum number v and we assume here P = B for more transparency. Similarly, the conductance 
can be found for d-wave symmetry, where the coherent gap is now a function of v, Ad, — Aocos20, </> is an angle 
along a constant energy contour. The DOS in the superconducting state is given by f8(| : 
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FIG. 5: (Colour online) Theoretical extrinsic conductance, Ea. (|16p for Ao — F, Ap — 2.7F and B = 2.65P. 
are indicated by the triangles and arrows respectively. Inset shows a representative STS spectrum of La2-i 
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Then, if the tail-width (F) is large compared with the coherent gap amplitude, F > Ao, Eq.(13) yields 



aNS oc A+ps{\eV\) [pAr(-eV) + PN{eV)] + 



1 arccos — r — B 1 — 

TT \Ao J \ Ao 



[p^{-eV)-pN{eV)], (16) 



where — 1 ± (^Q{~eV) + |) /P^. As aforementioned, the doping of impurities in the cuprates causes a 
bandtailing effect in the normal-state DOS. We have used a model DOS to reflect the shape of our DOS which is 
given by: 



PNjE) 
Pb 



1 + tanh ' ^ 



(17) 



This model pn{E) reflects the characteristic energy dependence of the DOS in disordered doped insulators, which is 
a constant ph above the two-dimensional band edge, and an exponent deep in the tail (see Fig|31). 

No matter what symmetry the superconducting gap is, the above equations capture all the unusual signatures of 
the extrinsic experimental tunnelling conductance in underdoped cuprates, such as the low energy coherent SG, the 
high energy PG and the asymmetry, see FiglS] In the case of atomically resolved STS one should replace the averaged 
Pn{E) in the above equations with a local bandtail DOS, pN^E^r), which is dependent on the position of the tip 
on the scanned area due to a nonuniform dopant distribution. As a result, the PG shows nanoscale inhomogeneity, 
while the low energy SG is spatially uniform as observed [i^ . Increasing doping level tends to diminish the bipolaron 
binding energy, Ap, since the pairing potential becomes weaker due to a partial screening of EPI with low frequency 
phonons [81.]. However, the coherent gap Ac, which is a product of the pairing potential and the squareroot of the 
carrier density [75| . can remain about a constant or even increase with doping, as observed [46l [6l|. 



IV. SS TUNNELLING 



For SS tunnelling there are different issues to address, in particular, tunnelling in mesas has indicated a nonzero 
conductance at zero voltage near and above Tc, also a negative excess resistance has been observed, along with the 
PG and SG. All of these features are accounted for in our theory. The Hamiltonian [t^ : 

vv' vv' V ^ / 

describes the tunnelling of a single polaron from one bosonic superconductor to the other. The first term has no 
involvement of bipolarons and describes the annihilation of a polaron in state v and creation of one in state v' . The 
second term involves the decay of a composed boson into two polarons, one remains in the same superconductor as the 
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FIG. 6: (Colour online) Cartoon demonstrating different single-particle tunnelling processes, the first is the annihilation of an 
electron in the superconductor on the left and the creation of a polaron on the right, the same as for NS tunnelling. This 
is described in the first part of the tunnelling Hamiltonian Ea. (|18p . The second illustrates a tunnelling process involving a 
bipolaron, where on the left a bipolaron is annihilated into two polarons, one of these moves into the polaron band on the left, 
the other tunnels to the superconductor on the right, as described by the second term of Eq.flS]). 



boson and is in state v (this is the time reversed state of v). The other tunnels into state v' in the other superconductor. 
The third term is the opposite to this with the annihilation of two polarons, one from each superconductor, that then 
combine to form a bipolaron. Again, only single-particle tunnelling is considered. Following the same procedure as 
for NS tunnelhng (applying the FDGR and Eq.dTUl)) yields: 



Iss{V) 



J2 + ^l^DU. - + eV- + ulvl,(U + - 1)] 

X [5{e^ + eV + e^,) - (5(e^ - eV + e^.)] 



+2ulvl, 



X [6ie^ -eV + e^')- S{e^ + eV + e^>)] 
/.(I - M] [H^. + eV -€,,)- 6(6, -eV- e,,)] • (19) 



The boson energy dispersion is again dropped here, assuming that bipolarons are sufficiently heavy for their bandwidth 
to be relatively small. 

First consider a clean bosonic superconductor, where there are no bandtails (F = 0), compare its I{V) with the SS 
I{V) in the BCS case for zero temperature when = 0. Refer back to Eg. dTTl) . the normal state DOS pn{(,) becomes 
a step function as F, which is the width of the bandtail, tends to zero, pn{0 ~^ ©(^ ~ ^p); Pn{^') — > B(^' — Ap). 
Thus, integrating with respect to ^ and ^' we have for the current: 



such that 



where /3 



Iss oc 2Q{eV - 2A) { (| + l) + |(eF - 2A) - [^{eV-A^-Al - ^A^ - 
{eVf 



2(A-Ae) 



eV + 2A 

e 



F(arcsin(l - ^), a) - {eV + 2Ac) [i^(arcsin(l - ^), a) - i;(arcsin(l - (3), a)] , 



(20) 



(21) 



'^'^ we take P = B, F{x, y) and E{x, y) are incomplete elliptic integrals of the first and 
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I{V) arb. units 




eV/A, 



FIG. 7: (Colour online) Our theory for SS tunnelling in a 'clean' bosonic superconductor, with no bandtail, compared to 
superconductor-superconductor tunnelling in the BCS theory [s^] at zero temperature for P = B, x/2 — 0.16, Ac = Abcs and 
A = 1.2Ac. 



second kind respectively. This is plotted in FigIT] against what is expected for SS tunnelling in the BCS theory [8! 

{eVf 



Iss ^e{eV-2ABcs) 



-K{a) - {eV + 2Abcs) [K{a) - E{a)] 



el/ + 2ABcs' 

where a = f^^fsf^f' ^l"^)' E{a) are complete elliptic integrals (see [s^ for more details.) 



(22) 



Finally consider the normal-state, where Aci/ = 0, which means €i, = \J Q + — \^A^ thus = 1 and = for 
positive ^jy, and = 0, = 1 for negative ^i,. We find: 



!-/.-/.')■ 



-eV + ) - Sic. + eV + e.' )] 



(23) 



Now, we can follow the same steps as Ref.[6l| and neglect temperature effects by approximating the distribution 
function f^, with a step function for temperatures near and above the transition temperature but sufficiently 

below the PG temperature T* ^ ^ > T > Tc, and for high enough voltages, \eV\ > ksT. Using the model 
normal-state DOS, Eq.dTT]) yields: 

-2 r a^(i + b^) 



IssiV) cx 



2(a2 ~ 1) 



In- 



1 



1 , a^ + b^^ 
^ In 



1 + 



B2(1-)-W2) , l + a%^ 
■ In 



In- 



P2(a2^4_i) a(l + 52) 2P2(a2_-64) a(l &2) 



(24) 



where a = e^r- and b = e~ . Near and above Tc, for high enough voltages eV > ksT, the conductance fits the 
experimental data and the gapped conductance is accounted for in underdoped mesas of Bi2212 50], as shown in FiglS] 
For these voltages the shape of the conductance is almost independent of the ratio B/P and x since the last term is 
much larger than all other terms in Eq.(24). The gapped conductance at smaller voltages can be accounted for by 
fully taking into account temperature effects in Ea. (l23|) . One can show that the quadratic term f„f„' in Eq. ((23| gives 
negligible contribution in the relevant temperature range, so differentiating the current with respect to the voltage 
yields the SS conductance as: 



ass oc 



sech 



e + - Aj 



sech 



C - - A, 



2P2 



PNi^O [fiO-fi-0] 



(25) 
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a (arb. units) 




FIG. 8: (Colour online) Approximate normal-state tunnelling conductance of bosonic superconductor (solid line (red)) with 
r = 3.2meV and Ap = 16meV 61] compared with experimental conductance 50| in mesas of Bi2212 (Tc = 95K at T = 85K). 



a (arb. units) 




10 20 30 40 50 60 70 

meV 



FIG. 9: (Colour online) Normal-state tunnelling conductance of a bosonic superconductor (solid lines) with fixed F = lOmeV 
and Bx^/P^ = 1.96 for each temperature, but changing Ap (see Fig llOl for temperature dependence of Ap) compared to 
experimental conductance [lOl for a few temperatures. 



This equation is plotted in Figl9] with F = lOmeV, the ratio B^x/2P^ is fixed at 1.96 and we change Ap to fit 
the experimental temperature dependence. Our theoretical results closely resemble those found experimentally by 
Krasnov [50| . Our comparison and other experiments 47, 59] suggests that the PG gradually decreases with increasing 
temperature, which could be explained by many-particle effects at significant doping. We suggest that with increasing 
temperature thermally excited mobile polarons screen the EPI, so the binding energy of the bipolarons drops [sij as 
with doping. 

Negative excess resistance below the transition temperature can be seen in cuprates [83| . Our theory can ac- 
count for this. Expanding Eas. ([T9)) and (|23p in powers of eV gives a zero bias conductance, for low temper- 
atures in the superconducting state this is 0-5(0) cx T^^ deps{e)^ cosh{e/2kBT)~^ , and in the normal-state 
crAr(O) cx T^^ d^p7v(C)^ cosh(^/2fcBT)~^. These integrals can be estimated to give respectively crs(O) oc 
T~-^ exp(— Ac/fcsT) for the s-wave coherent gap, or as{0) oc for the d-wave gap, and (TAr(O) oc exp(— T*/T). 
The latter expression is in excellent agreement with the temperature dependence of the mesa tunnelling conductance 
above Tc [sS (see also Ref.fsT']). Extrapolating this expression to temperatures below Tc yields the resistance ratio 
Rs/Rn oc Qy^-'I^B-T )/T (s_-^yave) or Rs/Rn oc exp{—T* /T)/T^ (d-wave). Hence in underdoped cuprates, where 
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FIG. 10: (Colour online) Temperature dependence of A reflects what is expected from experimental results [50 



T* > Ac/ks, the zero-bias tunnelhng resistance at temperatures below Tc is smaller than the normal state resistance 
extrapolated from above Tc to the same temperatures (i.e. the negative excess resistance), as observed (83| . 

V. SUMMARY 

Data from tunnelhng experiments is invaluable as it gives a huge insight into the low-energy excitations and thus 
the way high-temperature superconductors work. 

NS tunnelling in cuprates has indicated two energy scales, the first is the SG that vanishes above Tc- The second is 
the PG which remains a mystery, there is no general consensus as to what it is or why it exists. STM with cuprates 
has also shown that the tunnelling conductance of charge carriers in one direction (say tip to sample) is different 
to the tunnelling conduction in the other direction, this gives asymmetry in the tunnelling spectra. Again, there is 
no general consensus as to why this is the case. Remarkably, the position of the tip on the cuprate gives different 
tunnelling results. 

Intrinsic tunnelling and break junction experiments have also indicated a two-gap structure. Recent experiments 
have provided evidence that the PG is dependent on the temperature and decreases as the temperature increases 
above Tc- Another inexplicable feature of SS tunnelling has been the negative excess resistance, where the zero-bias 
tunnelling resistance that is extrapolated from above to below Tc is larger than that measured in the superconducting 
state below Tc- 

To have a successful theory, each of these puzzles should be accounted for in both NS and SS tunnelling. 

Our theory is based on the ab initio "LDA-I-GTB" band structures of charge transfer Mott-Hubbard insulators 
with the doping of impurities causing bandtails in the normal-state superconductor DOS. We make an extension to 
the BCS theory in the strong-coupling regime with bosonic (bipolaronic) carriers which are the real-space pairs of 
polarons. This theory has allowed us to describe the above unusual features of the cuprates. 

We greatly appreciate valuable discussions with Ivan Bozovic, Jozef Devreese, Kenjiro Gomes, Ruihua He, Viktor 
Kabanov, Vladimir Krasnov, Nikolai Kristoffel, Klim Kugel and Dragan Mihailovic. 
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